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ABSTRACT 

Measurements of pulsar flux densities are of great importance for understanding the 
pulsar emission mechanism and for predictions of pulsar survey yields and the pulsar 
population at large. Typically these flux densities are determined from phase-averaged 
“pulse profiles”, but this method has limited applicability at low frequencies because 
the observed pulses can easily be spread out by interstellar effects like scattering or 
dispersion, leading to a non-pulsed continuum component that is necessarily ignored 
in this type of analysis. In particular for the class of the millisecond pulsars (MSPs) 
at frequencies below 200 MHz, such interstellar effects can seriously compromise de¬ 
tectability and measured flux densities. In this paper we investigate MSP spectra 
based on a complementary approach, namely through investigation of archival con¬ 
tinuum imaging data. Even though these images lose sensitivity to pulsars since the 
on-pulse emission is averaged with off-pulse noise, they are insensitive to effects from 
scattering and provide a reliable way to determine the flux density and spectral indices 
of MSPs based on both pulsed and unpulsed components. Using the 74 MHz VLSSr 
as well as the 325 MHz WENSS and 1.4 GHz NVSS catalogues, we investigate the 
imaging flux densities of MSPs and evaluate the likelihood of spectral turn-overs in 
this population. We determine three new MSP spectral indices and identify six new 
MSPs with likely spectral turn-overs. 
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1 INTRODUCTION 

It is well established that there are two fundamental classes 
of pulsars: the normal (or “young”) pulsars, with spin pe¬ 
riods between roughly a tenth of a second and a few sec¬ 
onds and spin-down rates between 10~ 18 and 10~ 12 s/s; and 
the class of the millisecond pulsars (MSPs), with spin peri¬ 
ods below about 30 ms and spin-down rates below 10 -17 s/s 
(Lee et al. 2012). According to standard theory, this latter 
group is created through a “recycling” process where a nor¬ 
mal pulsar accretes mass from a main-sequence companion 
star through Roche-lobe overflow ( Alpar et al.||1982 l. This 
accretion process increases the spin of the pulsar but also 
affects a variety of emission properties, as summarised by 


Kramer et al. (19981. 


Studies of the spectral index of pulsars are of particu¬ 
lar importance to any attempts at understanding the pulsar 
emission mechanism; and to clarify any differences between 
MSPs and their slower counterparts. Even though only few 


MSP spectra were available, previous work (Kramer et al. 
1998 Toscano et al.|19'98 1 found no significant difference be¬ 


tween the spectral index distributions of MSPs and young 


pulsars, but Kramer et al. (1999) did find that MSP spectra 


tended to be somewhat less steep. Interestingly, recent sim¬ 


ulations by Bates et al. (20131 indicate that spectral indices 


available for slow pulsars are biased towards steeper spectra 
because most pulsars with determined spectral indices were 
discovered in early, low-frequency surveys, which are most 
sensitive to steep-spectrum objects. Arguably the same bias 
does not exist for the MSPs, which were mostly discovered 
in 1.4-GHz surveys and are less biased to steep-spectrum 
sources. Clearly, a full comparison of the populations will 
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require a significant increase in the number of spectral in¬ 
dices measured for MSPs. 

While spectral indices provide the simplest description 
of a spectrum, some pulsar spectra display more complex 
features. Of specific interest to our discussion are the spec¬ 
tral turn-overs that have been observed in some pulsars at 
frequencies below 400 MHz (e.g. Sieber|[l973 Kuzmin et al. 


|1978| |Izvekova et al.|1981|) and at frequencies near 1 GHz in 


a smaller sample ( Maron et al.|2000 l. 

Research into breaks in MSP spectra has been ham¬ 
pered by the small amount of data available; and incon¬ 
sistencies in this limited data set. |Erickson fe Mahoney| 
(1985) investigated the spectrum of the first known MSP 


(PSR J1939+2134 or B1937+21) and detected an indica¬ 
tion of a possible turn-over in the spectrum below 75 MHz. 


Subsequently, Foster et al. (1991) observed a sample of 


four MSPs, including PSR J1939+2134, at a variety of fre¬ 
quencies and found no significant deviations from power- 
law spectra in any of these pulsars. [McConnell et al.| 
(1996) published the spectrum of the brightest MSP known, 


PSR J0437—4715 and found clear evidence for a spectral 
break around 200 MHz. This pulsar remained the only MSP 
with a clear spectral turn-over until |Kuzmin fe Losovsky| 
([2001) published the largest sample so far, analysing the 
spectra of 30 MSPs at frequencies down to 100 MHz. Out of 
their large sample, only PSR J1012+5307 was found to have 
a turn-over in its spectrum. Most recently, observations of 
MSP J2145—0750 with the Long Wavelength Array in New 


Mexico (Ellingson et al. 2013) have shown very clear evi¬ 


dence of a flattening of the spectrum of this pulsar below a 
few hundred MHz ( Dowell et al.||2013 ), suggesting that the 
non-detection of a break by Kuzmin & Losovsky (2001) was 


caused by a lack of measurement points and sizeable uncer¬ 
tainties on their measurements. This brings the total num¬ 
ber of MSPs with suggested spectral breaks to four, out of a 
total sample of 33 MSPs with spectra down to ~ 100 MHz. 

As this historical overview shows, the spectra of MSPs 
have only been investigated down to ~ 100 MHz on a small 
number of sources; and when additional data were added to 
previously published spectra, the results were contradicting 
as often as they were confirming. This illustrates how lack 
of sensitivity; lack of observations; and the lack of access to 
a truly wide frequency range have hampered this Held. 

One effect that limits the accessible frequency range in 
which pulsar flux density measurements can be made, lies in 
the way most of these observations are conducted. Typical 
pulsar observations are performed using the phase-folding 
technique, in which the incoming data are averaged as a 
function of pulse phase, given the known rotational period of 
the pulsar. This technique is normally advantageous because 
it separates the pulsed signal from the off-pulse noise and 
thereby gains significant sensitivity. However, when investi¬ 
gating flux density measurements, this phase-folding tech¬ 
nique implicitly assumes that the off-pulse region (i.e. the 
phase-range where the pulsed emission reaches a minimum) 
only contains noise and no emission originating from the 


pulsar (Lorimer & Kramer 20051. This is not always true, 


because if a pulsar has only a small offset between its rota¬ 
tion axis and its magnetic axis, it is possible that the line of 
sight permanently crosses the polar emission region. In this 
case even the baseline level of the phase-folded profile would 


contain emission that should be attributed to the pulsar, 
but which gets ignored as it is unpulsed. 

At low frequencies this situation gets worse, since the 
pulsar emission is affected by a series of interstellar propaga¬ 
tion effects that tend to smear out the emission over pulse 


phase (see Rickett 1977 for a review). Because dispersion 


can nowadays be corrected for through coherent dedisper¬ 
sion, the most important of these ef fects is scattering , which 
has approximately a v ~ 39 scaling (B hat et al.| |2004) and is 
therefore particularly relevant at the low frequencies where 
the spectral turn-overs are observed to occur. 

An alternative approach to measuring spectra of MSPs 
is to determine the pulsar’s flux density from interferomet¬ 
ric imaging. Since these maps are not phase-resolved they 
do tend to have less sensitivity, but in this way the total flux 
density can be determinecQ (including any potentially un¬ 
pulsed component) and this flux density measurement is en¬ 
tirely insensitive to interstellar propagation delays like scat¬ 
tering. Moreover, several all-sky imaging surveys across a 
variety of observing frequencies have been released into the 
public domain, allowing spectrum investigations on a large 
sample of sources without any further requirement for ob¬ 
serving time. Such imaging surveys were already used by 


Kaplan et al. (1998) who investigated pulsar flux densities 


and scintillation based on the 1.4 GHz NVSS (NRAO VLA 
Sky Survey) and the 325 MHz WENSS (Westerbork North¬ 
ern Sky Survey), but they did not investigate these data 
for spectral breaks as they did not have access to a truly 
low-frequency survey. 


Since the analysis by Kaplan et al. (19981, the known 


sample of pulsars has more than doubled, and with the re¬ 
cent publication of the improved 74 MHz VLSSr (VLA Low- 
Frequency Sky Survey Redux, Lane et al. 20T2| , the sensitiv¬ 
ity of a similar investigation is significantly enhanced at the 
lower frequencies, providing significantly improved sensitiv¬ 
ity to spectral breaks over an unprecedentedly large number 
of sources. 

In this paper, we present the results from an investiga¬ 
tion into spectral turn-overs in MSPs, based on the 74 MHz 
VLSSr, complemented with the 325 MHz WENSS and the 
1.4 GHz NVSS. These three surveys are described in Sec¬ 
tion [2[ along with the selected source sample. Section [3] lists 
our detections and the obtained spectra; and in Section [4] 
we summarise our findings. 


2 IDENTIFICATION WITH APERTURE 
SYNTHESIS RADIO SURVEY IMAGES 

In order to investigate the spectra of MSPs, we investi¬ 
gated the data from the three imaging surveys discussed 
below, for detections of MSPs. In this context we defined 
MSPs as being pulsars with spin periods below 30 ms and 


1 Note that phase-resolved measurements of flux density typi¬ 
cally report the phase averaged flux density after subtracting a 
noise floor determined at off-pulse phases. This means that phase- 
resolved flux densities are directly comparable to interferometric 
flux density measurements. The only exception is the case where 
some part of the pulsar’s emission beam is continuously pointed 
towards us, as in this case only the pulsed fraction of the flux 
density is measured. 
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a time-derivative of the spin period P < 10~ 1 ' s/We ex¬ 
cluded pulsars in globular clusters to avoid obvious issues 
with confusion; and restricted our analysis to pulsars at de¬ 
clinations above —40°, as this is the lowest declination limit 
amongst the three surveys used. A list of candidate sources 
and th eir positions was taken from the ATNF pulsar cata¬ 
logue]^] (M anchester et al.|2005 ) and based on those positions 
the images from the three surveys were inspected for detec¬ 
tions. In some cases the positions from the catalogue are 
referred to epochs that are up to a decade or more after 
the data for the surveys were taken. This could cause the 
pulsar position to have changed due to the high velocity of 
pulsars. However, based on proper motion measurements of 
the MSPs in our sample (where available) or the average 
transverse velocity of the pulsar population (which is less 
than lOOkm/s, Hobbs et al.[2005), none of the MSPs in our 
sample was likely to have moved as much as a resolution 
element in any of the surveys used. 

Identification of our candidate sources was done 
through queries of the available source catalogue]/] and 
through visual inspection of the survey images, where a 3-er 
detection threshold coincident with the catalogue position of 
the pulsar was required to claim a detection. For the sources 
below ~ 5 <t significance, which were not in the online cata¬ 
logues, the position and flux density were subsequently de¬ 
termined through Gaussian fitting to the images. For all of 
our detections the offset between the fitted position and the 
most recent published position was well within the resolu¬ 
tion of the images. 

As mentioned, three imaging surveys were used. These 
were chosen based on their observing frequency, sensitivity 
and fraction of sky covered. The surveys used were: 

VLSSr: The original Very Large Array (VLA) Low- 
Frequency Sky Survey (VLSS) and its catalogue were re¬ 
leased in 2007 (Cohen et al. 20071. It covers 95% of the 


sky north of —30°declination at a frequency of 74 MHz with 
a resolution of approximately 80” and an RMS sensitivity 
of about lOOmJy/beam. The minimum integration time on 
each field is 75 minutes, consisting of three shorter obser¬ 
vations, each approximately 25 minutes, separated in time 
by at least one hour. Following improvements in various 
data reduction algorithms (particularly related to bright- 
source peeling, radio frequency interference removal and 
ionospheric calibration) a renewed release of the survey and 
images was made in 2012 (Lane et al.|2012 I. 

WENSS: The Westerbork Northern Sky Survey (WENSS), 
was performed with the Westerbork Synthesis Radio Tele¬ 
scope (WSRT), which consists of 14 parabolic antennas, be¬ 
tween 1991 and 1996 at a frequency of 325 MHz, covering 
the entire sky north of declination +28.5°with a limiting 1- 
(j flux density of approximately 3.6 mjy. This survey has an 


2 We also considered pulsars with as-yet undetermined spindown 
P, since these undetermined values are likely to be small. 

3 Available at http://www.atnf.csiro.au/research/pulsar/ 
psrcat 

The catalogue browsers for these surveys can be found at http: 
//www.cv.nrao.edu/vlss/VLSSlist.shtml, http://www.astron. 
nl/wow/testcode.php?survey=l and http://ww.cv.nrao.edu/ 
nvss/NVSSlist. shtml respectively. 


imaging resolution of 54" x 54"/ sin ( 5 ) with 5 the declination 
(Rengelink et al. 19971. 

NVSS: The NRAO VLA Sky Survey (NVSS) was per¬ 
formed with the VLA at a frequency of 1.4 GHz with the 
compact D and DnC configurations during 1993 and 1996, 
with some additional observations carried out in 1997. The 
NVSS is a continuum survey covering the entire sky north 
of —40°declination with a FWHM angular resolution of 45" 
and 1-a RMS flux density fluctuations of «0.45 mjy/beam 
( |Condon et al.| 1998]). 


3 RESULTS 

According to the ATNF catalogu]/] 130 sources satisfied the 
selection criteria laid out in Section [2] Out of these, we de¬ 
tected ten pulsars in the VLSSr data. A further six pulsars 
were expected to be detectable based on an extrapolation of 
the available higher-frequency data, but were not detected, 
providing strong evidence for a spectral break. The remain¬ 
ing 113 sources were not detected, but did not indicate a 
spectral turn-over. In some cases this was the case because 
their spectral index (as derived from higher-frequency data) 
is flat enough to make a non-detection at 74 MHz expected; 
while for other pulsars there are no flux density measure¬ 
ments at multiple frequencies available, making it impossible 
to derive a spectral index and therefore yielding no expected 
flux density at 74 MHz. The third and largest group (47 out 
of 113) were not detected because they lay outside the sur¬ 
vey area, because the VLSSr data for them were corrupted, 
or because their positions are not known precisely enough 
to allow a confident identification. 

In the following, we define the spectral index a as in 

S (y) oc v a , (1) 

with S (v) the flux density at frequency v. Since fitting non¬ 
linear functions in the presence of uncertainties is compli¬ 
cated and does not reliably quantify the often asymmetric 
uncertainties, we applied a Monte Carlo approach to the de¬ 
termination of our spectral indices. Practically this means 
we drew random realisations from the flux density measure¬ 
ments at all frequencies (using the published flux density val¬ 
ues and interpreting the measurement uncertainties as Gaus¬ 
sian and symmetric). A power-law function (Equation]!]) was 
then fitted to these simulated measurements through a lin¬ 
earised unweighted fit to determine realistic initial values 
which were used as the basis of a non-linear, weighted fit. 
The distribution of the slopes of these final fits (shown in 
figures |T] [5] and 0- then provided the most likely spectral 
index and its uncertainties. 

In the cases where a spectral turn-over was required, 
this was always clear without requiring a spectral fit, i.e. 
no borderline cases were found. In these cases a separate 
spectral slope was fitted to the lowest- and highest-frequency 
measurements, up to a break frequency that was visually 
identified. 

For the MSPs that were not detected in the VLSSr data, 
we used the same Monte Carlo approach on any available 
high-frequency data but extended the approach to provide 


5 Version 1.50 of the ATNF catalogue, used on 25 July 2014. 
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Table 2. Summary of spectral index measurements. This table 
contains the spectral indices derived from the measurements in 
Table [3 


Pulsar 

Spectral 

Turn-over 

a 

name 

index a 

frequency vto 

below i/to 

J0034—0534 

-2.64 ±0.05 

- 

- 

J0218+4232 

-2.41 ± 0.04 

- 

- 

J1810+1744 

9 rry+ 0.14 
z,o '- 0.16 

- 

- 

J1843—1113 

-3.08 ±0.10 

- 

- 

J1903+0327 

—2 01 +U ' U7 
z - ul - 0.05 

- 

- 

J1939+2134 

-2.59 ±0.04 

~ 74 MHz 

— 1 94 +0 06 
— 0.05 

J1959+2048 

9 70+U.O8 

z -' y -o.n 

- 


J2145—0750 

-2.60 ±0.04 

~ 400 MHz 

i qq+ 0.18 

i.oo—o i 7 

J2215+5135 

q 99+O.I5 
0.23 

- 

- 


a distribution of expected flux densities at an observing fre¬ 
quency of 74 MHz, assuming the power-law flux density scal¬ 
ing of Equation [l] This distribution was then compared to 
the 3 -<t upper limit on the flux density from the VLSSr im¬ 
ages to evaluate the likelihood of a spectral turn-over in 
these pulsars. 

In the following sections, we discuss our detections and 
the non-detections that may indicate spectral turn-overs. 
Specifically, in Section[3T]we discuss six detected MSPs that 
are consistent with previously published data; in Section |3.2| 
we present three detected MSPs that have their spectral 
index determined for the first time; and in Section |3.3| we 
discuss the six MSPs whose non-detections indicate likely 
spectral turn-overs. 


3.1 Six MSPs Detected and Consistent with 
Previous Publications 

In total, we detected 10 pulsars in the VLSSr. Four of these 
have no previous spectral index published and will therefore 
be discussed separately, in Section |3.2| The remaining six 
sources are all consistent with previously published spec¬ 
tral indices derived from higher-frequency data. Of these 
six, only PSR J0218+4232 was detected in the WENSS (the 
non-detections being primarily caused by the limited sky 
coverage of this survey) and the low luminosity of MSPs at 
higher frequencies made that only PSR J2145—0750 was de¬ 
tected in the NVSS. The list of the detected pulsars and all 
the flux density values used in our determination of spec¬ 
tral indices, is given in Table [l] the spectra of these pulsars 
are shown in Figure [T] and the derived spectral indices are 
given in Table [2] For a few sources a further discussion of 
our results is given below. 


3.1.1 PSR J0034-0534 

The spectral index of this pulsar has been measured twice 
before, by Toscano et al. (19981 who derived a spectral in¬ 
dex of a = —2.6 ± 1.0 and subsequently by |Kuzmin fe| 
Losovsky ( 2001|) who found a slightly shallower spectrum 


with a = —2.3 ±0.3. Our VLSSr detection lowers the lowest 
frequency in the spectrum from 111 MHz to 74 MHz and is 
in perfect agreement with the earlier Toscano et al. (19981 
measurements, but in significant disagreement with the flux 


density value of Kuzmin & Losovsky (2001). This disagree¬ 
ment indicates the Kuzmin & Losovsky (20011 flux density 
for this pulsar may have been underestimated. Excluding 
the Kuzmin & Losovsky (20011 value from our analysis, we 
derive a spectral index of a = —2.64±0.05. We note further¬ 
more that the folded pulse profile for this pulsar is very wide 
at these low frequencies]))] making it possible that some com¬ 
ponent of unpulsed, scattered flux density was unaccounted 
for in the Kuzmin & Losovsky (20011 measurement. 


3.1.2 PSR J0218+4232 


Navarro et al. (1995) discovered PSR J0218±4232 and re¬ 


ported a significant fraction of its radio emission was not 
pulsed, implying it is an aligned rotator and causing an off¬ 
set between its total flux density (as derived from imaging) 
and its pulsed flux density (as derived from phase-folded ob¬ 


servations). Stairs et al. (1999) subsequently determined the 


magnetic inclination angle to be consistent with 0°, based 
on a fit of the rotating vector model (RVM, see |Radhaki2 
ishnan & Cooke 19691 to observations at 410 and 610 MHz, 
thereby confirming the classification as an aligned rotator. 

We detected PSR J0218±4232 in the VLSSr data with 
a flux density of 2.9±0.4 Jy. In the WENSS, we detected the 
pulsar at 113mJy. We did not find the pulsar in the NVSS, 
which may be explained by scintillation (Rickett 1977j since 
the survey’s 1.35 mjy 3-cr sensitivity is close to the 1.5 mjy 
published flux density at 1.4 GHz. 

Combining our own measurements and those from 


Navarro et al. (19951 and Stairs et al. (1999) for both imag 


ing and phase-folded cases, we find no evidence for a spectral 
turn-over and determine a spectral index a = —2.41 ± 0.04 
for the total flux density, which is consistent with the spec¬ 
tral index for the phase-folded flux density: a = — 2.32±0.09, 
showing that the unpulsed flux density evolves consistently 
with the pulsed flux density (see Figure [I]. 


3.1.3 PSR J1939+2134 


This pulsar was originally known as 4C21.53 and was de¬ 
tected in many aperture synthesis images before being un¬ 
covered as the first MSP (then known as PSR B1937±21, 
Backer et al. 1982]. Consequently, its spectrum has been 
studied extensively, as can be seen from our summary in 
Figure |7] having flux density measurements from 4.8 GHz 
all the way down to 10 MHz. An initial controversy about 
the presence of a spectral turn-over could be explained by 
the very slight nature of the turn-over and the very low fre¬ 
quencies at which it occurs. 

We detected PSR J1939±2134 in the VLSSr data with 


a flux density of 17.94 ± 2.15 Jy. This measurement is con¬ 
sistent with previously published flux density measurements 
and confirms the slight spectral turn-over first indicated by 


Erickson & Mahoney (1985). The WENSS did not cover 


this position and a confident detection in the NVSS was 
made impossible by a coincident extended source of emis¬ 
sion slightly north of the pulsar position. 


6 See, e.g. the European Pulsar Network Database at 
http://www.jb.man.ac.uk/research/pulsar/Resources/epn/ 
browser.html 
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Table 1. Summary of flux density measurements used. This table contains the flux densities used for our determination of spectral 
indices (summarised in Table |2j, combining measurements from the VLSSr, WENSS and NVSS surveys with those already available 
from literature. In cases where the measurement uncertainty was not clear from the publication, we assumed a 30% uncertainty. These 
measurements are indicated by an asterisk. 


J1810+1744 

J1816+4510 
J1843—1113 

J1903+0327 


J1939+2134 

(B1937+21) 


J1959+2048 

(B1957+20) 


Pulsar name 

frequency 

(MHz) 

flux density 
(mjy) 

Reference and Notes 



J0034—0534 

74 

1620 ± 210 

This work. 






103 

250 ± 120 

Kuzmin & Lo 

sovsky 


2001 

), significantly lower flux density, not used in fit 


436 

17 ± 5 

Toscano et al. 

(1998 




660 

5.5 ±0.5 

Toscano et al. 

(1998 

) 



1400 

0.61 ± 0.09 

Toscano et al. 

(1998 

) 



1660 

0.56 ±0.09 

Toscano et al. 

(1998 



J0218±4232 

102 

270 ± 150 

Kuzmin &; Lo 

sovsky 


2001 

- 

(pulsed) 

410 

24 ± 7* 

Stairs et al. 

L999 





410 

20 ± 10 

Navarro et al. 

iff 

95 




606 

8 ± 5 

Navarro et al. 

(1955 




610 

15 ± 7 

Kramer et al. 

(1998 




610 

11 ± 3* 

Stairs et al. (] 

999 





1410 

0.9 ±0.2 

Navarro et al. 

(19 

J5 | 




1410 

0.9 ±0.2 

Kramer et al. 

(19- 

18 



J0218±4232 

34.5 

(4± 1) x 10 4 

Dwarakanath 

6£ Udaya Shankar 

(1990 

- 


74 
151 
325 
325 
608 
1400 
74 
350 
74 
74 
1400 
74 
1400 
2000 
5000 
74 

various 

700 

1400 

2695 

3000 

4850 

8350 

74 

318 

430 

606 

1490 

1400 


2940 ± 360 
660 ± 330 
113 ±4 
150 ± 50 
26 ± 7.8 

1.5 ±0.5 
1090 ± 150 

20 ± 4 
690 ± 110 
820 ± 194 
0.10 ±0.02 
440 ± 100 
1.3 ±0.4 
0.62 ±0.05 
0.09 ±0.02 
17940 ± 2150 
various 
63 ± 19 
13 ±5 
2.0 ±0.4 

1.6 ±0.7 

1.0 ±0.2 

0.039 ±0.015 
1880 ± 240 
60 ±40 
26 ± 13 
8 ± 5 
0.4 ±0.2 
0.4 ±0.2 


Navarro et al. 
This work. 


(1995 i 


Navarro et al. 

(1995 

Navarro et al. 

(1995 

Navarro et al. 

(1995 


Hessels et al. 
This work. 

This work. 


12011 ) 


Hobbs et al. 
This work. 


120041 


Champion et al. 

(2008 

Champion et al. 

(2008 

Champion et al. 

(2008 


See 

Erickson & Mahoney i 

1985) 

Manchester et al. 

(2013 


Manchester et al. 

(2013 


Kramer et al. (1999) 


Manchester et al. ( 2013 i 

Kramer et al. (1999) 


Kowalinska et al. (2012) 


Fruchter et al. 

(1990 

) (retrieved from their Figure 2) 

Kruchter et al. 

(1990 

(retrieved from their Figure 2) 

Fruchter et al. 

(1990 

(retrieved from their Figure 2) 

Fruchter et al. 

(1990 

(retrieved from their Figure 2) 

Kramer et al. ( 

1998:1 



3.1-4 PSR J1959+2048 


PSR J1959±2048 (also known as PSR B1957±20) is the first 
known eclipsing pulsar, discovered at 430 MHz by [Fruchter| 


et al. (1988). The eclipses occur when the pulsar passes be¬ 


hind its degenerate, 0.1 -M@ companion, which is being ab¬ 
lated by the pulsar wind. At the discovery frequency, the 
pulsar is eclipsed ~10% of the time, but this eclipse fraction 
scales strongly with wavelength ( Stappers et al.|20 01). 

We detected PSR J1959±2048 in the VLSSr data with 


a flux density of 1.88 ± 0.24 Jy. This area was not covered in 
the WENSS and we did not detect the pulsar in the NVSS 


either, possibly due to eclipses. Our measurement extends 
the spectrum further down, indicating the absence of spec¬ 
tral breaks. 


3.1.5 PSR J2145-0750 


Dowell et al. (2013) recently observed PSR J2145—0750 with 


the Long Wavelength Array (LWA) in the frequency range 
41 to 81 MHz. Combining their measurement with archival 
higher-frequency data, they fitted the pulsar’s spectrum to 
the following, somewhat unconventional, functional form: 


© 2014 RAS, MNRAS 000, |T|{To] 
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Table 1 - continued 


Pulsar name 

frequency 

(MHz) 

flux density 

(mjy) 

Reference and Notes 


J2145—0750 

various 

various 

See 

Dowell et al. (2 

013 



74 

354 ± 91 

This work. 






102 

480 ±120 

Kuzmin & Losovsk 

y(2001 



102.5 

90 ± 50 

Maloleev et al 

. 120C 

10 ) 



410 

46 ± 14* 

Stairs et al. 

i 

999 




430 

50 ± 15* 

Bailes et al. 


994) 




436 

100 ± 30 

Toscano et al. 

(199 

81 



610 

19 ±6* 

Stairs et al. 

i 

999) 




660 

36 ±4 

Toscano et a . 

(1998) 



700 

16 ±9 

Manchester et 

al. ( 

2013 

) 


1400 

7± 0.9 

Toscano et al. 

(199 

4 



1400 

2.6 ±0.4 

This work; significantly lower flux density, not used in fit 


1414 

6.6 ±2.0* 

Stairs et al. 

(jl999) 




1510 

8 ± 2 

Kramer et al. 

(199S 

il 



1520 

10 ±3* 

Bailes et al. 

< 

994.) 




1660 

5.4 ± 0.8 

Toscano et al. 

(199 

81 



2695 

2.2 ±0.5 

Kramer et al. 

(1991 

1 ) 



3100 

1.4 ±0.5 

Manchester et al. ( 

2013 

) 


4850 

0.4 ± 0.1 

Kramer et al. 

(1999 

> 



4850 

0.44 ± 0.03 

Kijak et al. 

1997 




8350 

0.080 ± 0.037 

Kowalinska 

et al. (1 

2012 

) 

J2215+5135 

74 

800 ±120 

This work. 






350 

5 ± 1.5* 

Hessels et al. 

2011 

) 






J1959+2048 



10 2 10 ; 
Frequency [MHz] 


1(T icr 

Frequency [MHz] 



Figure 1. Spectra for the six detected pulsars with existing spectral information. Top row: spectra for three MSPs with no sign of a 
spectral break. Bottom row: spectra for the three MSPs with more complex spectra. The spectra for PSRs J1939+2134 and J2145—0750 
show spectral turn-overs; and the spectrum for PSR J0218+4232 is consistent but at a higher level in imaging data than in phase-folded 
data as a consequence of the unpulsed emission (see Section |3.1.2| for details). The insets show the probability density functions for the 
spectral index distribution for each pulsar; in these insets the most likely spectral index and its 1 -a confidence interval are indicated by 
vertical bars. For the three pulsars on the bottom row, two different spectra were fitted. These two fits were based on the flux density 
measurements shown with dashed symbols; and flux density measurements shown in full symbols. In practice, for PSR J0218+4232, this 
means that dashed measurements are from phase-folded work, while the full markers indicate flux densities derived from imaging data. 
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Spectral Turn-Overs in Millisecond Pulsars 7 


SM = s„(i) 


1 / \ mi+m2 log(i//z/ r ) 


( 2 ) 


with m i the spectral index at the lowest frequencies, m 2 the 
spectral curvature and v r the rollover frequency, which was 


relatively arbitrarily chosen by Dowell et al. (20131 to be 
about 730 MHz. 

We detected PSR J2145—0750 in the VLSSr image with 
a flux density of 354 ± 91 mjy. This flux density measure¬ 
ment is consistent with the measurement by |Dowell et al.| 
|2013|, who determined a flux density of 390±65±35 mjy at 
73 MHz, where the given uncertainties are for measurement 
noise and systematic uncertainties, respectively. We also de¬ 
tected this pulsar in the NVSS image, having a flux density 
of 2.6 mjy, which is significantly less than the 7.0±0.9 mjy of 


Toscano et al. (19981. This discrepancy could be explained 


by scintillation, as indicated by the long-term monitoring 


of Manchester et al. (20131, who found the flux density of 


PSR J2145—0750 at 1.4 GHz to have an RMS scatter of no 
less than 12.5 mjy. The pulsar’s position was not included 
intheWENSS. 

Because no flux density measurements are published 
for this pulsar between ~100MHz and ~400MHz, the ex¬ 
act break frequency is ill-determined (if there is indeed a 
particular break frequency rather than the continuous turn¬ 
over advocated by [Dowell et al. 2013). Nevertheless, it is 
clear that at frequencies below 100 MHz the spectrum is 
flatter than above 400 MHz. Based on our own and the 
archival data listed above, we determine spectral indices of 
— 1.33^017 and —2.60± 0.04 respectively. 


3.2 Three New Spectral Indices 


Four pulsars we detected in the VLSSr only had previous 
flux density measurements at a single frequency (or none at 
all), implying we were able to derive three spectral indices 
for the first time. These three spectra are shown in Figure[2] 
and the flux density measurements and spectral information 
are summarised in tables [l] and [2] respectively. Interestingly, 
two of the four pulsars in this class, PSRs J1810+1744 and 
J2215+5135 were found in a follow-up survey of unidenti¬ 
fied 7 -ray sources in the Fermi all-sky map (Hessels et al. 


2011) and three out of four (PSRs J1810+1744, J1816+4510 


and J2215+5135) inhabit eclipsing binary systems (Hes- 
|sels et aL[|2011| |Kaplan et al. ||2012[ ) . We cannot investigate 
possible spectral turn-overs for these pulsars, as we have 
less than three flux density measurements per pulsar. For 
PSR J1816+4510 specifically, we do not determine a spec¬ 
tral index either, as the only flux density measurements cur¬ 
rently available are derived from uncalibrated survey obser¬ 
vations and are therefore unreliable (see [Stovall et al. 2014[ 
and their Figure 6 in particular^ 

Finally, given their detections at low frequencies, these 
pulsars are by design biased towards steep spectral indices. 
The measured spectral indices are between a = —2.57 and 
a = —3.22, confirming that these pulsars constitute a partic¬ 
ularly steep-spectrum part of the pulsar population: of the 


7 If we were to use the |Stovall et al.| |2014| flux density values, 
an extreme spectral index of —3.76 would result, well in excess 
of the steepest pulsar spectral index currently published (-3.5 for 


PSR J0711+0931, 


Lommen et al 


. 2000 ). 



-3 -2-10 1 

Spectral index 


Figure 3. Histogram of all known pulsar spectral indices. The 
main figure shows all known spectral indices of pulsars, the inset 
shows the same plot for the class of the millisecond pulsars only. 
Our three new measurements are shown in black. 


328 spectral indices available in the ATNF pulsar catalogue, 
only 9% have equally steep spectra as these three MSPs (see 
Figure [3|. 


3.3 Six MSPs with a High Probability of Spectral 
Turn-Over 

Based on fluxes at higher frequencies (either available from 
literature, or derived from the NVSS or WENSS data), a 
spectral index can be derived and used to predict the flux 
at 74 MHz. For six undetected pulsars this predicted flux 
lay above the 3-cr detection limit of the VLSSr. This may 
in some cases be explained by scintillation or because the 
higher-frequency flux density estimates could be unreliable 
('Levin et al. 2013). With the exception of those caveats, a 
non-detection can also imply a spectral turn-over. For these 
unexpectedly non-detected pulsars, we determined the prob¬ 
ability of a spectral turn-over using the Monte Carlo ap¬ 
proach outlined at the very start of this section. Table [3] 
provides the list of these pulsars, along with the spectral 
indices derived from the published flux densities and the 
likelihood of spectral turn-over following from our analysis. 
For PSR J1744—1134 we also used our NVSS detection at 
3.5 ± 0.5 mjy. The spectra are shown in Figure [4] 


4 CONCLUSIONS 

There are two fundamentally different ways to determine 
pulsar flux densities. These are through observations aver¬ 
aged modulo the rotational period of the pulsar (“phase- 
folded” measurements); and through interferometric imag¬ 
ing observations. The former type has higher sensitivity be¬ 
cause of the separation of pulsed signal and unpulsed noise, 
but the latter approach is unaffected by scattering in the 
interstellar medium (which especially affects low-frequency 
data) and does provide a correct measurement of the to¬ 
tal pulsar flux density even in the case of aligned rotators, 
where only a part of the emission is pulsed. 


© 2014 RAS, MNRAS 000, [T]{To] 
































































8 M. Kuniyoshi et. al. 




Figure 2. Three new MSP spectra obtained from our own imaging flux densities and published data. The inset panels show the likelihood 
distribution of the spectral indices, with the vertical lines indicating the most likely value and its l-a interval. 


Table 3. Summary table for seven MSPs with a high probability of having a spectral break. Given are the pulsar name, the spectral 
index derived from previously published flux densities, the probability of turn-over and the literature references for the flux density 
measurements. Note these probabilities do not account for scintillation or underestimated uncertainties in the published measurements. 
References are as follows: (1) Kuzmin &; Losovskyj (|2001 1 ; (2) |Kramer et al.~| (|1998[ > ; (3) |Kramer et aI7| (|1999 1; (4)|Lorimer et ak|il995~| l; 
(5) |Toscano et al.| (| 1 998ty ; (6) Manchester et a.1. | (|2013[ l ; (7)|Stairs et al.| (|1999[ t ; (8) |Lorimer et al.| ( | 1996| ) ; (9) |Deneva et al.| l |2012 l; (10) 
|Boriakoff et al ,| ( | 1 983^ ; (11) |Lommen et al.| ( |2000[ ) ; (12) |Kaplan et al.| ( | 1998[ ) ■ 


Pulsar 

Spectral 

Turn-over 

References 

name 

index a 

likelihood 


J0030+0451 

—2 14+ 0 - 25 
Z ' 14t -0.52 

77% 

i, ii 

J1640+2224 

o 1 q+O. 12 

97% 

1, 2, 3 

J1643—1224 

oii +0.03 
Z ‘ 11 -0.04 

100% 

2, 3, 4, 5, 6, 12 

J1744-1134 

1 qc+0.06 

—0.08 

100% 

This work and 1, 2, 3, 5, 6, 7 

J1911—1114 

-2 45+ 0 - 06 
z -^°-0.05 

100% 

1, 2, 5, 7, 8 

J1949+3106 

q i q+ 0.24 
C *- i -0.17 

100% 

This work and 9 

J1955+2908 

—2 29 +0,34 

0.25 

82% 

2 and 10 


In this paper, we have used archival imaging data to 
determine the spectra of MSPs. We identified pulsars in 
the 1.4 GHz NVSS, the 325 MHz WENSS and, most impor¬ 
tantly, the 74 MHz VLSSr. The VLSSr is of particular im¬ 
portance because it extends the investigated spectra to the 
lower frequencies, where phase-folded observations are often 
hampered by scattering effects and where investigation of 
potential spectral turn-overs becomes possible. 


We only obtained a few detections in the NVSS, mostly 
because of the limited sensitivity of this survey compared to 
the low flux densities of MSPs at 1.4 GHz. In the WENSS 
our detection rate was lowered by the survey’s limited sky 
coverage. The VLSSr did yield ten new pulsar flux density 
measurements; and for six sources a non-detection in the 
VLSSr provided significant evidence for spectral turn-overs 
at lower frequencies. Such evidence is hard to obtain from 
phase-folded observations since those cannot distinguish be¬ 
tween pulsars that are too weak to be detected and pulsars 
that are bright but heavily scattered. 


In total, we extended the spectra of 16 pulsars to fre¬ 
quencies below 100 MHz and found evidence for turn-overs 
in eight of these. Five of the eight pulsars with spectral turn¬ 
overs were previously reported to not have a spectral turn¬ 
over, based on higher-frequency data of |Kuzmin fe Losovsky| 
(2001), suggesting in these cases the turn-over occurs near or 
just below 100 MHz. Following our work, a total of 39 MSPs 
have had their spectra investigated at frequencies down to or 


below 100 MHz; and 10 out of those show evidence for turn¬ 
overs - a much higher fraction than previously thought. 

The largest analysis of spectra for young pulsars to 
date, was published by Maron et al. (2000). They found 
there were two separate populations of pulsars with spec¬ 
tral breaks. About 10% of the population displayed spectral 
steepening at frequencies around or above 1.4 GHz, whereas 
a far smaller number of pulsars (on the order of a percent) 
showed spectral turn-overs at frequencies around 100 MHz. 
The frequencies of the spectral turn-overs found in our sam¬ 
ple are comparable to the latter, low-frequency turn-over 
sub-population. The fraction of the population displaying 
turn-overs (about one in four), however, is similar to the 
fraction of young pulsars showing high-frequency turn-overs 
(~ 10%). This lends some credence to the postulate that 
MSP spectra are fundamentally identical to the spectra of 
young pulsars, albeit shifted to lower frequencies, as first 
proposed by (Kramer et al. 1999]). 

Future complements to our work will be provided by a 
host of new SKA pathfinder arrays that operate at the low¬ 
est frequencies continuously observable from Earth. All three 
of these pathfinders (LOFAR, LWA, MWA) are now capa¬ 


ble of performing phase-resolved pulsar observations ( Stap- 
pers et al. 2011 Dowell et al. 2013| Bhat et al. 2014 1 and 


will therefore carry spectral investigations to lower frequen¬ 
cies, even for fainter pulsars that are undetectable in most 
imaging surveys. Furthermore, all-sky surveys by these three 
telescopes, combined with phase-resolved observations may 
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Figure 4. Spectra for six MSPs with a high probability of having a spectral turn-over. The arrows indicate the 3-cr upper limits in the 
74 MHz VLSSr maps. The insets in the top-right of each sub-figure show the probability distribution of the spectral index based on the 
higher-frequency data referenced in Table [3] while the inset in the lower-left of each sub-figure shows the probability density function 
of the 74 MHz flux density based on this distribution of spectral indices, along with the detection probability at 74 MHz assuming the 
spectrum does not have a spectral break. The vertical lines in the top-right plot show the most likely spectral index value along with its 
1 -ct uncertainty interval; the vertical line in the lower-left plot shows the 3 -a detection limit of the VLSSr (approximately 300 mjy). 


shed some light on the prevalence and strength of interstel¬ 
lar scattering in the sample of MSPs that are undetectable 
in phase-resolved observations. 
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